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Abstract. We study acoustic white holes in a steadily flowing atomic Bose-Einstein 
condensate. A white hole configuration is obtained when the flow velocity goes from 
a super-sonic value in the upstream region to a sub-sonic one in the downstream 
region. The scattering of phonon wavepackets on a white hole horizon is numerically 
studied in terms of the Gross-Pitaevskii equation of mean-field theory: dynamical 
stability of the acoustic white hole is found, as well as a signature of a nonlinear 
back-action of the incident phonon wavepacket onto the horizon. The correlation 
pattern of density fluctuations is numerically studied by means of the truncated- Wigner 
method which includes quantum fluctuations. Signatures of the white hole radiation 
of correlated phonon pairs by the horizon are characterized; analogies and differences 
with Hawking radiation from acoustic black holes are discussed. In particular, a short 
wavelength feature is identified in the density correlation function, whose amplitude 
steadily grows in time since the formation of the horizon. The numerical observations 
are quantitatively interpreted by means of an analytical Bogoliubov theory of quantum 
fluctuations for a white hole configuration within the step-like horizon approximation. 
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1. Introduction 

The recent experimental advances in the creation and manipulation of atomic Bose- 
Einstein condensates are suggesting these systems as ideal candidates where to 
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experimentally address fundamental questions in quantum hydrodynamics [1]. In 
particular, a lot of attention has been devoted to configurations showing regions of 
super-sonic flow, where the dynamics is the richest and quantum effects should be 
easiest to observe. Among the most remarkable recent observations, we may mention 
the emission of a wake of Bogoliubov phonons by a moving defect via an analog of 
the Cerenkov effect [2], the shedding of solitons in a one-dimensional configuration [3], 
the onset of dynamical instabilities in the periodic potential of optical lattices [4]. In 
the meanwhile, theoretical investigations have addressed general questions about the 
stability of superflow [5] and have identified the mechanisms underlying the shedding of 
solitons [6] and vortices [7]. Analytical aspects of the one-dimensional flow have been 
also investigated in detail [8] . Most of this physics is accurately described within a mean- 
field approach, where the dynamics of the Bose-condensed atomic cloud is described in 
terms of the Gross-Pitaevskii equation for the condensate wavefunction [1, 9]. 

In the wake of the formal analogy [10, 11] between the propagation of quantum 
fields on a curved space-time background [12] and the propagation of sound on 
inhomogeneously moving fluids, a strong interest is presently being devoted to the 
quantum features of the condensate hydrodynamics. A most intriguing prediction 
concerns the emission of correlated pairs of phonons by the horizon of an acoustic 
black hole configuration via a mechanism which is a condensed-matter analog of the 
well-celebrated Hawking radiation from gravitational black holes [13]. Following the 
proposal of [14] , the key signature of this effect was identified in the correlation function 
of density fluctuations: the correlation between the Hawking phonon and the partner 
falling into the black hole results in a peculiar, tongue-shaped feature in the correlation 
pattern. This prediction was confirmed in an ah initio numerical simulation of the 
dynamics of an atomic BEC [15]. A first experimental attempt to create an acoustic 
black hole in an atomic condensate was reported in [16]. 

As we schematically show in Fig.l, an acoustic white hole configuration is obtained 
by simply reversing the direction of flow of a black hole one: atoms go from an upstream, 
inner region where the flow is super-sonic to a downstream, external region where 
the flow is sub-sonic: within a naive hydrodynamic picture, dragging by the moving 
condensate forbids low-energy phonons from crossing the horizon and penetrating the 
inner region of the white hole. Questions about the stability of such acoustic white 
hole configurations has lately attracted a lot of interest, and different answers have been 
proposed to this problem, ranging from instability [17], to stability [18], to intermediate 
behaviours depending on boundary conditions [20]. Issues related to the quantum 
emission of phonons by configurations involving white hole horizons have been addressed 
in [18, 19, 21, 22]. 

The present paper reports a comprehensive theoretical investigation of the physical 
properties of acoustic white holes in atomic Bose condensates. Focussing our attention 
on a specific, yet realistic configuration which is most suitable to numerical and 
analytical study, we address the problem of the dynamical stability of white hole 
configurations and then the physical properties of the white hole radiation emitted by 
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the horizon. In order to be able to safely isolate the dynamics of the white hole horizon 
from spurious effects due to the finite size and/or the multiply connected geometry of the 
system [23], all numerical calculations have been performed using absorbing boundary 
conditions for the Bogoliubov phonons so to model the most relevant situation of an 
infinitely long condensate. 

In Sec. 2 we introduce the physical system and we discuss the main features of the 
dispersion of Bogoliubov excitations in respectively the super- and sub-sonic regions: 
this will be the essential ingredient to physically understand the numerical observations 
that are presented in the following of the paper. Among the different possible choices of 
the density and flow velocity pattern that result in a white hole horizon, we concentrate 
our attention on a configuration that appears most robust against spurious effects. 

The dynamical stability of the white hole configuration is assessed in Sec. 3, where 
the scattering of a phonon wavepacket on the white hole horizon is studied by means 
of the Gross-Pitaevskii equation of mean-field theory. In addition to the reflected 
wavepacket in the outward direction, a pair of short wavelength wavepackets penetrate 
inside the white hole: as a consequence of the strong blue shift experienced while 
approaching the horizon, low energy hydrodynamic modes are in fact transformed 
into high-momentum single particle excitations that are able to enter the white hole. 
One of the two wavepackets consists of negative-norm Bogoliubov modes and results 
from the same mixing of positive and negative norm modes that underlies Hawking 
radiation. Dynamical stability of the configuration is demonstrated by the fact that 
no exponentially growing deformation of the horizon develops under the effect of the 
incident wavepacket \. The only visible back-action effect consists of the Bogoliubov- 
Cerenkov emission of zero-frequency phonons [2] in the upstream direction by the slight 
deformation of the horizon; a specific discussion of this feature is given in Sec. 3. 2. As 
a function of time, the amplitude of this pattern eventually tends to a finite amplitude 
value, roughly proportional to the square of the incident wavepacket amplitude. 

The effect of quantum fluctuations is studied in Sec. 4 using the same truncated 
Wigner method as in [15]. In analogy to what was done for acoustic black holes, our 
key observable is the correlation function of density fluctuations. Several characteristic 
features are identified and related to the correlated emission of phonon pairs by the 
white hole horizon into different pairs of modes. In spite of some similarities with the 
Hawking emission from black holes, the physical properties of the emission turn out 
to be significantly different: for this reason, we have chosen to refer to it under the 
different name of white hole emission. While Hawking radiation mostly consists of low- 
k excitations, the phonons that are emitted into the white hole have a large momentum 
comparable to the healing length of the condensate; furthemore, the checkerboard 
pattern that appears in the correlation function for points located inside the white 
hole keeps growing in time with a characteristic logarithmic or linear law depending on 

| A similar procedure for more complex black hole/white hole configurations has confirmed the black 
hole laser phenomenon [21, 24], i.e. the onset of an exponentially growing perturbation of the density 
profile. This will be the subject of a forthcoming publication. 
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the initial temperature. This is to be contrasted to the late-time stationary state of 
Hawking radiation from black holes. 

An analytical understanding of the white hole radiation is developed in Sec. 5 based 
on the Bogoliubov theory of dilute Bose Einstein condensates within the step-like horizon 
approximation of [25]. Dynamical stability of the white hole configuration is assessed by 
checking that complex frequency eigenmodes of the Bogoliubov-de Gennes equations are 
indeed absent. A theory of quantum fluctuations is developed in terms of the S'-matrix 
describing phonon scattering on the white hole horizon. This provides an analytical 
expression for the spectral distribution of the white hole emission: in contrast to the 
thermal character of Hawking radiation from black holes, the phonons that are emitted 
outside the white hole have a flat spectral distribution. The same approach is then used 
to evaluate the different contributions to the correlation pattern: the resulting analytical 
expressions reproduce in a fairly accurate way all the features that are observed in 
the numerical calculations. The steady growth in time of the checkerboard pattern is 
related to an infra-red divergence of the density correlation function at late times: the 
logarithmic/linear growth is recovered by imposing a suitable infra-red cut-off to the 
integral. Conclusions are finally drawn in Sec. 6. 

2. The physical system and the theoretical model 
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Figure 1. Left panel: sketch of the acoustic white hole configuration. Right panel: 
sketch of an acoustic black hole configuration. The two configurations are related by 
a reversal of the condensate flow speed vq- 



A sketch of the physical system we are considering is shown in the left panel of 
Fig.l. An elongated atomic Bose-Einstein condensate is steadily flowing along an atomic 
waveguide in the x direction. A suitable modulation of the confinement potential and/or 
of the atom-atom scattering length creates a white hole horizon at x = separating a 
upstream x < region of super-sonic flow from a downstream x > one of sub-sonic 
flow. In the gravitational analogy [10, 26, 27], the upstream region corresponds to the 
interior of the white hole, while the downstream region corresponds to the external space. 
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For comparison, a sketch of the acoustic black hole configuration obtained reversing the 
flow speed is shown in the right panel §. 

To simplify the theoretical description, the transverse confinement is assumed to 
be tight enough for the transverse degrees of freedom to be frozen [1] and the system 
dynamics to be accurately described by a one-dimensional model based on the following 
second-quantized Hamiltonian, 

h 2 



% — I dx 



V¥(x) Vtf(ar) + V(x, t) ¥(x) V(x) + 

+ ¥{x) ¥{x) %{x) V(xj\ . (1) 

Here, ^f(x) and W(x) are atomic field operators satisfying Bose commutation rules 
^(x')] = 8(x — x'), m is the atomic mass, V(x) the external potential, and g(x) 
is the effective one-dimensional atom-atom interaction constant [1]. At mean-field level, 
the condensate dynamics is described by the one-dimensional Gross-Pitaevskii equation 
(GPE) [1] 

ih Tt = S + t] * + 9{X > t] ^ * (2) 

for the macroscopic condensate wavefunction ip. In Sec. 3, this equation will be used to 
study the propagation of classical perturbations on top of the white hole configurations: 
full inclusion of the nonlinearity will allow us to study the dynamical stability of 
the configuration beyond the linearized Bogoliubov theory. In the following Sec. 4, 
the so-called truncated Wigner method will be adopted to go beyond the mean- 
field approximation underlying (2) and include the effect of quantum and thermal 
fluctuations: these are described in terms of stochastic noise on the initial value of 
the condensate wavefunction, which then evolves according to the same GPE (2). 
Expectation values of observables are obtained as averages over the stochastic noise. 
We refer to [15] for a detailed description of the numerical technique. 

For the sake of simplicity, we shall focus our attention onto a specific configuration 
which is most suitable for both analytic and numeric investigation. Initially, the 
condensate is assumed to have a spatially uniform density no and a spatially uniform 
flow speed v along the positive x direction. The external potential and the (repulsive) 
atom-atom interaction constant are also uniform and equal to respectively V(x) = Vd 
and g(x) — g^ > 0. Around t = to, a step-like spatial modulation is applied to both the 
potential and the interaction constant by suitably modifying the transverse confinement 
potential and/or the atom-atom scattering length via an external magnetic field tuned 
in the vicinity of a so-called Feshbach resonance [1]: within a short time a t , V and g in 
the upstream x < region are brought to their final values V u and g u , while their values 
in the x > region are kept equal to the initial ones Vd and g^- The transition region 
around x = has a characteristic thickness a x ; away from it, both V and g quickly tend 
to their asymptotic values V u ^ and g u d . 

§ Note that this configuration differs from the one considered in [15, 25] by a simple spatial inversion 
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To avoid the development of spurious instabilities due to the reflection of Bogoliubov 
excitations at the edges of the integration box, suitable absorbing boundary conditions 
are implemented for the Bogoliubov modes. In practice, this is done by including a 
restoring force on the wavefunction that damps out Bogoliubov excitations and brings 
it back to a pure condensate. For this additional term not to interfere with the 
physics under investigation, it is spatially restricted to small regions at the edges of 
the integration box which are sufficiently far from the white hole horizon. 

In order to reduce the impact of competing processes such as back-scattering of 
condensate atoms and soliton shedding from the potential step [6], the external potential 
V is chosen to exactly compensate the spatial jump in the Hartree interaction energy 
fht,d = 9u,d n o, i-e. 

V d + n d = V u + n u . (3) 

In this way, the plane wave 

ip(x,t) = y/no exp[i(kox — uot)} (4) 

with hk /m = vq and Uq = hk^jlm is for all times a solution of the GPE (2). 

The white hole configuration that is the subject of the present paper is characterized 
by the chain inequality c u < Vq < c d relating the speed of sound in respectively 
the upstream and downstream asymptotic regions c u ^ = a/ l^ u ,d/ m and the (spatially 
uniform) flow speed vq. The dispersion of the Bogoliubov modes describing in the 
laboratory frame the propagation of weak disturbances on top of the flowing condensate 
has the usual form 



^2/^2^2 \ 

u(k) = n u ,d( k ) + v o k = ± \ ^ ( "2^- + 2mc l,d J + v o k (5) 

in respectively the upstream and downstream regions and is illustrated in the upper (a,b) 
panels of Fig. 2. For future use, we have defined Q Ujd (k) as the frequencies measured in 
the upstream and down stream co moving frames at rest with the condensate. 

The main feature of the white hole configuration consists of the fact that long 
wavelength phonons are forbidden from penetrating into the inner region of the white 
hole. As a consequence of dragging by the super-sonic flow, the group velocity turns 
out to be positive Vq ± c u > in the laboratory frame for both k > and k < 
phonons. This behaviour has to be compared to the black hole case shown in the lower 
panels of Fig.2: in this case, the direction of the condensate flow is reversed and long 
wavelength phonons cannot escape from the inner, super-sonic region. It is however 
crucial to note that the above black and white hole behaviours are restricted to \ow-k 
Bogoliubov excitations of hydrodynamic nature. Indeed, large-A; Bogoliubov excitations 
have a single particle character and propagate at a faster group velocity. As a result, 
they are able to propagate in both directions, that is to penetrate into a white hole and 

|| Note that this convention differs from the one that is usually adopted in relativity, where the modes 
are labelled according to their left- or right- moving character with respect to the light/sound cone, i.e. 
their propagation direction in the frame comoving with the condensate [18, 19]. 
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Figure 2. Dispersion of Bogoliubov excitations in the asymptotic regions far away 
from the horizon. The labels of the modes indicate the up- or down-stream region with 
respect to the condensate flow and the in- or out-going direction of their propagation 
direction (as predicted by their group velocity) with respect to the horizon ||. Both the 
frequency ui(k) and the group velocity of the modes are measured in the laboratory 
frame. Panels (a,b): white hole configuration with a rightward flow v$ > 0, as sketched 
in Fig. 1(a). White hole parameters: v /c u = 1.5, v Q /cd = 0.75. Panels (c,d): black 
hole configuration with a leftward flow vq < 0, as sketched in Fig. 1(b). Black hole 
parameters: |«o|/ c d = 1-5, |i>o|/c« = 0.75. 



escape from a black hole. The transition from the hydro dynamic to the single-particle 
regime occurs at a characteristic wavevector given by the inverse of the healing length, 

t,u,d = h/ mc Uid . 

In terms of the gravitational analogy, the point separating the upstream region of 
super-sonic c u < t>o flow from the downstream one of sub-sonic v o < Q flow then acts as 
the horizon. The k < hydrodynamic phonons that propagate in the sub-sonic region 
in the leftward direction towards the horizon get slowed down as they approach the 
horizon and can never penetrate into the supersonic region. As a result, they pile up 
in the vicinity of the horizon. Analogously, the group velocity of k < hydrodynamic 
phonons in the super-sonic region tends to zero as the horizon is approached. As a 
result, they get also piled up on the horizon. The mechanism of this pile-up effect is 
graphically illustrated in the upper-left panel of Fig.3. As propagation on a stationary 
background conserves the phonon frequency u in the laboratory frame, the pile-up effect 
is in both cases accompanied by an increase of the wavevector k and a corresponding 
blue shift of the phonon frequency Q(k) in the comoving frame. While this blue-shift 
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Figure 3. Sketch of the mode propagation around a horizon of the white (left panels) 
and black (right panels) hole type. Upper panels refer to the case of linear dispersion. 
Lower panels include the super-luminal contribution to the dispersion. 



effect continues for ever in the case of dispersionless (relativistic) fields giving rise to 
excitations at arbitrarily high k values, in a Bose-Einstein condensate is restricted to 
wavevectors smaller than the inverse healing length by the peculiar super-luminal shape 
of the Bogoliubov dispersion shown in Fig. 2. 

As the white hole configuration is related to the black hole one by a simple time 
reversal, the corresponding Bogoliubov dispersions shown in the upper (a,b) and lower 
(c,d) panels of Fig. 2 are connected by a k — > —k symmetry. From the figure, it is 
immediate to see that under this operation the sign of the group velocity is reversed 
and the in- and out-going character of the modes are exchanged. Another, related 
difference involves the zero- mode defined by the u(kz) = condition. Such mode only 
exists in the super-sonic region and propagates opposite to the condensate flow: while in 
a black hole it propagates towards the horizon, in a white hole it propagates away from 
the horizon into the inner region. As a result, a perturbation of the white hole horizon 
from the Hartree condition (3) may result into the Landau-Cerenkov [2] emission of 
phonons into this zero-energy mode of wavevector 

2m 



kz = — f~\jvl-cl (6) 
comparable to the inverse healing length and group velocity 

vz = (7) 
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In the following of the paper, we shall see how this fact plays a crucial role in the physics 
of the white hole configuration. 

3. Dynamical stability of a white hole configuration: mean-field theory 

3.1. Scattering of phonon wavepackets on the horizon 



to 



IXlQ- 
5xlO" J 



-5x10 

-lxlO' J 
lxl 0' J 

5x1 0" J 



K> -5x10 



CO 



1 1 1 1 1 


1 1 ] 1 1 1 1 1 

d-in 


1 1 . 


- Kd , 1 . 1 , 


; —*t$^Pfi»— 




1 



: 4tik|i 


1 1 1 1 1 1 


— 11, t = 1250 


1,1,1," 




Figure 4. Scattering of a phonon wavepacket on a white hole horizon. The different 
panels show snapshots of the density modulation 8n(x) = n{x) — n Q at subsequent 
times. The insets provide a closer look at some most interesting features. White hole 
parameters as in Fig.2(a,b) with a x /£d — 0.5. Wavepacket parameters k wp ^ d = 0.2, 
o'wp/id — 80. The qualitative behavior remains unchanged if smoother horizons are 
considered with a x /£d ^ 1- 



In the previous section we have seen that a condensate wavefunction in the plane 
wave form (4) is at all times a solution of the GPE (2). Still, one has to assess the 
dynamical stability of this solution. A simple way to answer to this question consists 
of solving the Bogoliubov-de Gennes equations for the spectrum of linear perturbations 
on top of the flowing condensate and looking for complex frequency modes. Pioneering 
calculations in this direction were performed in [20] for different black and white hole 
configurations and for different boundary conditions. A detailed investigation of black 
hole laser instabilities in combined black and white hole configurations was reported 
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in [21, 22]. Application of this method to the white hole case is briefly discussed in 
Sec.5.1. 

In the present Section, we shall follow a different method which consists in studying 
the wavefunction evolution under the full GPE (2) for an initial condition slightly 
perturbed from the plane wave solution (4). This approach takes into account the 
nonlinear couplings that are neglected in the Bogoliubov theory and that are responsible 
for interesting back-reaction effects. In all calculations, the white hole is assumed to 
be already formed well before the experiment is performed so that the condensate 
background can be safely considered as a stationary one. 

In particular, we consider the simple case where a phonon wavepacket modulation 
has been imprinted on top of the homogeneous solution (4). The wavepacket has a 
carrier at k wp and a Gaussian envelope of width £ wp much larger than the wavelength, 
k wp £ W p ^ 1- Initially, the wavepacket is centered at a position x° far from the horizon, 
\ x wp\ ^ ^w P - In wavevector space, the phonon wavepacket has a narrow Gaussian shape 
around k wp of width Ak wp = l/£ wp <C k wp . The amplitude of the wavepacket is chosen to 
be weak enough to be in the linear regime; this condition has been checked by verifying 
that the amplitude of the reflected and transmitted wavepackets scales linearly with the 
amplitude of the incident wavepacket. A remarkable feature stemming from nonlinear 
effects beyond this approximation is discussed in Sec. 3. 2. 

Snapshots of the wavepacket evolution at subsequent times are given in Fig. 4 for 
the most illustrative case in which the wavepacket is initially at x° wp > outside the 
white hole and propagates towards the horizon on the d m mode with a wavevector 
k wp < well within the hydrodynamic region, \k wp \ <C Q 1 . Its frequency distribution 
is concentrated around the carrier frequency u wp ~ (v — q) \k wp \. As the wavevector 
distribution is concentrated in the hydrodynamic region k^d «C 1 where the dispersion 
is almost linear, the wavepacket gets weakly distorted while propagating towards the 
horizon. Dynamical stability of the white hole configuration is confirmed by the fact 
that no exponentially growing perturbation of the horizon is triggered by the arrival of 
the incident wavepacket. 

Within a strict hydrodynamic picture where the high momentum m™' modes were 
completely neglected, one would expect that no excitation can propagate across the 
horizon into the white hole as no leftward propagating modes are available inside the 
white hole at u wp . As a consequence, the incident wavepacket either accumulates at 
the horizon or is reflected. This picture is oversimplified as the combination of the 
blue shift at the white hole horizon and the super-luminal character of the Bogoliubov 
dispersion may lead to a significant transmission into modes with large wavevectors 
in the vicinity of rkkz- In [18], it was predicted that the amplitude of these transmitted 
wavepackets is much larger than the one of the reflected one even for very smooth white 
hole horizons. 

These expectations are confirmed in our numerical calculations shown in Fig. 4. The 
reflected wavepacket is actually small, but clearly visible in the insets of the two latest 
time panels. As expected, it propagates away from the horizon at the speed Cd + vq. 
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In agreement with the condition that the carrier frequency oo wp has to be conserved 
after scattering, the carrier wavevector of the reflected wave packet is smaller than the 
incident one by a factor (v — Cd)/(vo + q). 

At relatively short times after impinging on the horizon (see e.g. the fx^t = 2000 
panel of Fig. 4), there is a single transmitted wavepacket that propagates into the white 
hole with a group velocity close to vz- As one can see in the left inset, its shape is 
characterized by a carrier at a wavevector kz and a slower modulation. This latter can 
be interpreted as the result of interference between the u° ut and u% ut components of 
wavevector approximately A; u o Ut ~ k z + ou wp /vz and k u out ~ — kz + u wp /vz, see Fig. 2(a). 
From the Bogoliubov expression for the density perturbation [9], one immediately 
recognizes oscillations at ±k u out and ±A; n out. The slow modulation corresponds to the 
beating of the ±/c M out and ^fk^m components, and has a wavevector q^eat ~ 2u wp /\vz\ ~ 
2(cd — vq) \k wp \/\vz\- The quantitative agreement of this Bogoliubov prediction for qb eat 
with the period of the slow modulation extracted from the [i^t = 2000 panel of Fig. 4 
confirms our interpretation. 

At later times (e.g. /i^t = 3500), the curvature of the Bogoliubov dispersion makes 
the u° ut and u™ 1 wavepackets to eventually separate in space. As a result, the slow 
modulation due to their interference disappears and one is left with two Gaussian 
wavepackets. The presence of the transmitted wavepacket on the negative norm u™ 1 
out-going mode can be interpreted as the result of a stimulated Hawking emission [25] . 
This classical counterpart of Hawking radiation was recently investigated in experiments 
using gravity waves in water tanks [28, 29]. 

Had we considered an incident wave packet with a smaller wave vector k wp , the 
time needed for the two transmitted wave packets to separate would be longer as the 
group velocities associated with the two roots u° ut and u™ 1 would be closer. In the 
\k W p\ limit, the two wave vectors k°f" and k°^ tend to ±k z and the wavepackets 
never separate. 

For the sake of completeness, it is interesting to compare these results to the case 
of a black hole configuration. As before, we consider a Gaussian wavepacket centered at 
x wp > outside the black hole which propagates towards the horizon on the u m mode 
with a wavevector k wp < well within the hydrodynamic region, \k wp \ <C As one 
can see in the snapshots shown in Fig. 5, the u m incident wave packet scatters on the 
horizon and splits into a reflected wavepacket on the u out mode and a pair of transmitted 
ones on the modes labelled as d° ut and d?, wt in Fig.2(c,d): while the stronger and faster 
wavepacket corresponds to the transmitted positive-norm dl ut mode, the weaker one 
corresponds to the negative-norm d™ 1 mode and results from the same mixing of positive 
and negative modes that underlies Hawking radiation [25]. In contrast to the white hole 
case, all out-going wave packets have significantly different group velocities and separate 
in space almost immediately; furthermore, the carrier wavevectors of all wavepackets 
tend to zero in the \k wp \ — > limit. The reason of these remarkable differences is easily 
understood in terms of the dispersion curves shown in Fig.2. 
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Figure 5. Scattering of a phonon wavepacket on a black hole horizon. The different 
panels show snapshots of the density modulation 8n(x) = n(x) — uq at subsequent 
times. Black hole parameters as in Fig.2(c,d) with a x /^ u = 0.5. Wavepacket 
parameters k wp ^ u = 0.05, a wp /^ u = 80. 



3.2. Back-reaction of the wave packet on the horizon 

The previous discussion of phonon scattering on the white hole horizon has highlighted 
a quite similar physics as compared to the black hole case: an incident wave packet 
splits in several reflected and transmitted wavepackets whose carrier wavevector is fixed 
by energy conservation arguments. As we are considering weak incident wavepackets, 
the amplitude of the reflected and transmitted scales linearly with the incident one. 
However, in addition to these expected wave packets, a striking unexpected feature 
is visible in the latest time snapshot of Fig.4 right inside the white hole horizon: a 
stationary density modulation with a spatially flat envelope and a fast oscillating carrier. 
At late times, the amplitude of this density modulation saturates to a finite value. A 
similar effect was very recently reported for a classical white hole analog using gravity 
waves on the surface of a water tank [28, 29]. 

Its physical origin can be understood in terms of the zero-mode in the Bogoliubov 
dispersion: the vanishing w-frequency of this mode explains the stationarity of the 
pattern and the observed wavelength of the modulation coincides with the value kz of 
the zero-mode wavevector in the upstream region. Still, the microscopic rectification 
process that is responsible for the excitation of a zero-frequency mode from a finite 
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Figure 6. Log-log plot of the amplitude of the zero-mode modulation as a function 
of the peak incident wavepacket amplitude. The zero-mode modulation amplitude 
is extracted from late-time snapshots of simulations of the same configuration as 
considered in Fig. 4. The black dotted line is a fit with a square law 8n z = C Sn^, p . 

frequency wavepacket requires further investigation as it is not allowed within a 
linearized Bogoliubov theory on a stationary background. 

To assess its origin in terms of some nonlinear effect, we have repeated the 
same numerical simulation of the GPE with different amplitudes of the incident 
wavepacket. The results are summarized in Fig.6: the late-time amplitude of the zero- 
mode modulation scales indeed proportionally to the square of the incident wavepacket 
amplitude. A simplest physical interpretation can be put forward as follows: the incident 
wavepacket slightly distorts the wavefunction in the horizon region breaking the Hartree 
condition (3). The resulting perturbation of the Hartree potential is then responsible 
for the continuous emission of phonons via the Landau- Cerenkov emission discussed 
in [2, 6] for a uniform super-sonic flow hitting a defect. 

This interpretation is confirmed by the numerical observation that no such feature 
was ever observed for black hole configurations: in this latter case, the super- sonic region 
is in fact located downstream of the defect so that the zero-mode propagates towards 
the horizon. This fact is clearly visible in Fig. 2(c). As a result, the zero-mode can not 
be excited by a perturbation localized in the horizon region of a black hole. 

To summarize, our numerical simulations of the Gross-Pitaevskii equation show 
that the white hole configuration is dynamically stable at mean-field level even beyond 
the linearized regime considered so far in the literature. Even when it is excited by a 
sizable incident wavepacket, the condensate does not develop any exponentially growing 
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perturbation: the amplitude of the scattered wavepackets is roughly proportional to the 
incident amplitude and the amplitude of the zero-mode perturbation tends at late times 
to a stationary value proportional to the square of the incident wavepacket amplitude. 



4. The effect of quantum fluctuations 

After having assessed in the previous Section the dynamical stability of a white 
hole configuration at the level of mean-field theory, we now proceed to investigate 
the consequences of quantum fluctuations, and in particular the radiation that is 
continuously emitted by the white hole horizon once this has been formed. Our 
investigation will make use of the same truncated- Wigner numerical technique used 
in [15]. We refer to this paper for a more detailed discussion about technical details. 

At the initial time t — 0, the condensate is assumed to be spatially homogeneous at 
density no and to be flowing at a constant speed vq; both the interaction constant and 
the external potential are initially flat and equal to g& and Vd, respectively. The white 
hole horizon is formed at the time to by ramping g and V in the upstream x < region 
to their final values g u and V u . This ramp takes place within a time a t . Eventually, g 
and V tend to a arctan-shaped spatial profile of characteristic thickness o~ x . 

Quantum and thermal fluctuations on top of the initial homogeneous condensate 
can be treated by means of the standard Bogoliubov theory of dilute Bose gas. Within 
the Wigner method, this translates [30] into a stochastic initial condition with a suitable 
Gaussian noise added on top of the plane-wave initial wavefunction (4) of mean-field 
theory. The variance of the Gaussian noise is determined by the initial temperature T . 
Exception made for Fig. 9(b), all numerical data shown in the figures refer to the case of a 
vanishing initial temperature T = 0. The time evolution of the stochastic wavefunction 
starting from its initial state follows the deterministic GPE evolution (2) including the 
chosen form of the nonlinear interaction constant g(x,t) and of the external potential 
V(x,t). Expectation values of symmetrically ordered observables at any later time are 
obtained by taking the corresponding stochastic average over the ensemble of evolved 
wavefunctions. 



4-1- Signatures of white hole emission in the correlation pattern of density fluctuations 

In [15] it was shown that a good deal of information of the physical properties of Hawking 
radiation can be extracted from the correlation pattern of density fluctuations at equal 
times, defined as 

(&(x) &(x')if(x')if(x)) 



G^(x,x') = 



(W(x)V(x)) (W(x')V(x')) 
(n(x)n(x')) 1 



. S(x — x'). (8) 
(n(x)) (n(x')) (n(x)) K J K ' 

In this section, we follow the same line for the case of a white hole, trying to isolate the 
peculiar features that characterize the white hole counterpart of the Hawking radiation 
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Figure 7. Color plots of the rescaled density correlation function (no£d) x [G^ 2 ' (x, x') — 
1] at two successive times njt — 90 (a) and 160 (b) after the switch-on of the white hole 
horizon. The (i), (m), (iv) labels identify the main features discussed in the text. 
The solid magenta and red lines indicate the directions along which the cuts shown in 
Fig. 8 (a, b) are taken. The white rectangle indicate the region where the checkerboard 
amplitude shown in Fig. 9 is measured. The blue dashed line in panel (b) indicates 
the analytical position of the axis of feature (iii) . The horizon is formed within a time 
Hd &t — 10 around /idto = 50 and has a spatial thickness (J x /£,d = 0.5. White hole 
parameters: vq/c u = 1.5, vo/cd = 0.75. The initial temperature is To = 0. 
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Figure 8. Cuts of the rescaled density correlation function at times fx\t — 160 (dotted 
lines), 180 (solid lines) taken along respectively the (a) magenta (x' — 63.5^) and (b) 
red (x — x' = 15.5^) lines of Fig. 7. The dashed line in (a) indicates the prediction 
(26) for the central position of the d — u lt2 correlation discussed in Sec. 5. 3. 



from black holes. Looking at the Bogoliubov dispersions shown in Fig.2(a,b), we expect 
from energy conservation arguments that correlated pairs of quanta can be emitted from 
the horizon either into w° ut and modes or into the d ont and modes: the negative 
energy of the u™ 1 partner can in fact compensate the positive energy of the u° ut or 
d ont ones, so to give a zero total energy of the pair. As it happened in black holes, 
these processes correspond to clearly distinct features in the correlation pattern. An 
analytical understanding of this white hole radiation will be provided in Sec. 5. Color 
plots of a suitably normalized density correlation G^(x, x') are shown in Fig. 7. For 
numerical convenience, all density correlation plots have been smoothened out with a 
Gaussian spatial filter of size l a vj^d — 1- The main features that are visible in the plots 
can be classified as follows: 

(i) The strong and negative correlation strip on the x = x' main diagonal results from 
the usual many-body antibunching due to repulsive interactions and has no relation 
with the presence of an horizon. 

(ii) A system of stripes parallel to the main diagonal appears in the x < region inside 
the white hole as soon as the horizon is formed. As time goes on, the fringes move 





Figure 9. Time evolution of the peak-to-peak amplitude of the checkerboard pattern 
measured within the region indicated by the white rectangle in Fig. 7. The upper (a) 
panel shows a semilog plot for a zero initial temperature To — 0. The lower (b) panel 
shows a linear plot for different temperatures fcsTo//irf = (black), 0.1 (red) and 0.2 
(blue). The cyan straight lines are guides to the eye which emphasize that the growth 
is approximatively logarithmic (linear) in the zero (finite) temperature cases. 



away from the main diagonal at an approximately constant speed and eventually 
disappear from the region of sight. We have checked that their presence and shape 
does not depend on the presence of the horizon, but rather on the speed of the 
time-modulation: the slower the ramp (i.e. the longer at), the weaker the fringe 
amplitude. On the basis of [31], these fringes can be interpreted as a consequence 
of the dynamical Casimir emission of pairs of counterpropagating phonons at all 
points where the interaction constant is modulated in time, g = ga — > g u - It is 
interesting to note that in the white hole case the dynamical Casimir phonons are 
able to travel back to the horizon. As a direct consequence of this fact, feature (ii) 
extends also into the x > 0, x' < and x < 0, x' > quadrants. As the scattering 
of dynamical Casimir phonons incident off the horizon may interfere with the white 
hole emission, it is then important to ensure that the intensity of the dynamical 
Casimir emission is efficiently suppressed by a suitably long formation time a t . 

(in) Spatially oscillating correlations appear between points right outside and right 
inside the horizon. The correlation pattern has a sinusoidal dependence as 
a function of the position of the point located inside the horizon, while the 
dependence on the outside point is much smoother. The wavevector of the fringes 



Acoustic white holes in atomic BECs 



18 



is approximately equal to the zero-mode wavevector kz introduced in the previous 
section. As time goes on, these fringes extend into a larger area further away from 
the horizon, but their amplitude at each given point eventually tends to a finite 
value, constant in time. This last fact is better visible in the cuts taken along the 
magenta line that are shown in Fig.8(a). Remarkably, the envelope of the fringe 
system shows a slow modulation, with a minimum approximately located along the 
blue dashed line of Fig. 7(b), x/v z = x' /(v + c d ), which suggests that this feature is 
related to the emission of pairs of quanta by the horizon into the ui^ and d modes. 
An analytical discussion of this feature will be provided in Sec. 5. 3. 

(iv) A checkerboard pattern appears in the correlations between pair of points located 
inside the white hole. As time goes on, the checkerboard pattern extends further 
away from the horizon. The spatial wavevector of the checkerboard is approximately 
equal to kz in both x, x' directions, which suggests that it is due to the emission 
of pairs of quanta by the horizon into the u\$ modes. In contrast to feature (Hi) 
and to what was always observed for black holes in [15], its amplitude at any given 
point does not tend to a stationary value, but keeps on steadily increasing. This 
remarkable fact is visible in the cuts of the correlation pattern shown in Fig. 8(b) 
and, even more clearly, in the plot of the checkerboard amplitude as a function of 
time shown in Fig.9(a,b). For a vanishing initial temperature T = 0, the growth 
appears to follow a logarithmic law. For a finite initial temperature T > 0, the 
growth appears to be instead linear with a slope proportional to T . An analytical 
interpretation of all these features will be discussed in Sec. 5. 4. 

4-2. Comparison with Hawking radiation from black holes 

For the sake of completeness it can be useful to briefly review the main properties of 
the Hawking emission from a black hole [15, 25]. This should facilitate the reader in 
appreciating the analogies and differences with the white hole radiation. As before, 
this study is based on the correlation pattern of density fluctuations, some examples of 
which are shown in Fig. 10. Several features can be identified: 

(i) As expected, the strong and negative correlation strip on the x — x' main diagonal 

is almost identical in the white hole and black hole configurations. 
(ii) The moving pattern of fringes due to the dynamical Casimir emission of phonons by 
the time-modulated interaction constant is also very similar in the black and white 
hole configuration. The main difference is purely geometrical: as a result of the 
different relative position of the super- and sub-sonic regions, the low- A; dynamical 
Casimir phonons are not able to travel back to the horizon and are fully dragged 
into the black hole. 

(Hi) A symmetric pair of negative correlation tongues extend from the horizon point for 
points x and x' located on opposite sides from the horizon. Their axis is located 
along the x/(vq + q) = x' / (vq + c u ) straight line. While their length linearly grows 
with time since the horizon formation time, their maximum height remains almost 
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Figure 10. Color plots of the rescaled density correlation function (no£d) x 
[G^ (x, x') — 1] at two successive times [i u t — 40 (a) and 160 (b) after the switch-on of a 
black hole horizon. The (i), (ii), (Hi), (iv) labels identify the main features discussed 
in the text. The horizon is formed within a time \i u a t = 0.5 around fi u to — 2.5 
and has a spatial thickness er x /£ u ~ 0.5. The calculation has been performed using 
the truncated- Wigner method. Same system parameters as in Fig.2(c,d). The initial 
temperature is T = 0. 



constant in time. Their physical origin lies in the stationary emission of correlated 
pairs of quanta by the horizon into the modes labelled as u out and c?2 Ut i n Fig-2(c,d): 
while the d<i partner falls into the black hole, the u partner flies away giving rise to 
the standard Hawking emission. 

(iv) Another pair of positive correlation tongues extends from the horizon point for both 
x, x' < inside the horizon. Their axis is located along the x/(vo + Cd) = x'/ (v — Cd) 
line, which suggests an interpretation in terms of the stationary emission of 
correlated pairs of d° ut and d^ quanta by the horizon. 

(v) Another pair of weaker tongues appears in the x < 0, x' > region along the 
x/(vo — Cd) = x'/(vo + c u ) line, which suggests an interpretation in terms of 
correlations between the d° nt and u out modes. 

A general fact of all (iii)-(v) features is that none of them shows any significant 
modulation: this is due to the fact that the wavevector of all out-going phonon modes 
tends to zero in a black hole configuration for low frequencies. This is to be contrasted 
to the white hole case where the wavevector of modes tends to ±kz which is of the 
order of the inverse healing length in the supersonic region. 

To summarize, in this Section we have described the properties of the quantum 
emission of phonons from the white hole horizon. As in the black hole case, the peculiar 
quantum correlations between the consituents of each pairs result into long distance 
correlations of the density fluctuations. Several novel features have been highlighted: 
as the quantum emission into the white hole consists of short wavelength excitations 
comparable to the healing length, the correlation pattern shows a fast sinusoidal 
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modulation. Furthermore, the amplitude of the checkerboard pattern for pairs of points 
located inside the white hole keeps on growing in time. This is in stark contrast with 
the stationary character of Hawking radiation from black holes. The consequence of this 
steady growth of density fluctuations for the stability of the white hole configuration 
requires a detailed study of the back-action effect of quantum fluctuations onto the flow 
and will be the subject of further work. 



5. Analytical theory 

In this Section we sketch an analytical framework offering a physical interpretation of 
the numerical observations discussed in the previous Sections. More details will be given 
in future publications, in particular in Carlos Mayoral's PhD thesis [32]. Our approach 
is based on the step-like horizon model first introduced at the level of mean-field theory 
in [20] and then developed into a full theory of quantum fluctuations and Hawking 
radiation in [25]. Here, we will show that this theoretical model is able to explain most 
of the numerically observed physics and to provide quantitatively accurate predictions 
for most relevant observables. 



5.1. Dynamical stability of white holes 

In this first subsection we address the issue of the dynamical stability of the white hole 
configuration. The numerical observations reported in Sec. 3 indicate that the chosen 
configuration is dynamically stable. Here we confirm this conclusion by analytical 
means. Our analytical approach is inspired by [20] and gets to similar conclusions. 
For analytical simplicity, we restrict our attention to the same configuration used in 
the numerical calculations, that is a BEC in the plane wave state (4) with a spatially 
constant density n and flow speed v ] the step-like white hole horizon is created by 
means of a sudden jump at x = of both the nonlinear interaction constant g(x) and 
the external potential V(x). 

Dynamical stability of this configuration against perturbations of the horizon region 
is assessed by looking for eigenmodes of the Bogoliubov-de Gennes equations [1, 9] with 
complex frequency: the system is unstable if at least an eigenmode a exists with a 
frequency u a such that Im[u) a ] > 0. As we are interested in instabilities that start in 
the horizon region, boundary conditions have to imposed that the mode wavefunction of 
the perturbation decays to zero at infinity on both sides of the horizon. At the horizon, 
the mode wavefunction has to be continuous along with its first spatial derivative. 

More precisely, we are looking for eigenmodes of the Bogoliubov-de Gennes 
equations: 

h 2 d 2 

hcuu(x) = - ^g^2 U ( x ) + V(x)u(x) + 2g(x)\(f) \ 2 u(x) + <f>l(x)v(x) - fJ,u($) 
ft 2 d 2 

huv(x) = - — -—^v(x) — V(x)v(x) — 2g(x)\(f) \ 2 v(x) — (f) (x)u(x) + fxv(x)(10) 

Zi III (J X 
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Within the spatially homogeneous upstream x < region, the eigenmodes can be 
expanded into plane waves as follows 

u (x) = a ^u (11) 

i 

0*0® e*°", (12) 

i 

where the sum over i runs over the Bogoliubov modes at frequency Q, whose complex 
wavevector ku satisfies the (complex) dispersion relation 

h\Q - k u v ) 2 = ^ + 2g u \4> \ 2 ) (13) 

2m \ 2m I 

describing Bogoliubov modes in a moving fluid at vq. As usual, the amplitudes u u and 
v u are determined by inserting the plane wave ansatz (11-12) into the Bogoliubov-de 
Gennes equations (9-10). Differently from the standard cases [1, 9], the complex value 
of ui makes the u, v coefficients to be complex as well. 




Re[tfl] 

Figure 11. Log-scale color plot of the determinant whose zeros define the 
asymptotically bound modes of a white hole configuration with a step-like horizon. 
X- and y-axis are the real and imaginary parts of the complex frequency G). As the 
determinant has no zeros in complex u) plane, no asymptotically bound modes exist. 
This proves a fortiori that the white hole is dynamically stable. Same white hole 
parameters as in the previous figures. 

The boundary condition that the mode wavefunction has to decay to zero for 
x — » — oo restricts the sum to those Bogoliubov modes whose wavevector Imffci''] < 0: 
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out of the four solutions of (13), only two can then be retained %. Analogous conditions 
hold for the downstream, x > region. The matching condition at the horizon imposes 
that both u(x) and v(x) have to be continuous, as well as their first spatial derivatives. 

The boundary conditions at infinity plus the matching conditions at the horizon 
reduce our problem to the detemination of four amplitude coefficients (two modes per 
u, d side) that have to satisfy four matching conditions. Non-trivial solutions (the so- 
called asympotically bound modes (ABM) [21, 22]) exist only for special values of the 
complex frequency Co that make the determinant of the corresponding linear system to 
vanish. As one can see in Fig. 11, this never happens for the white hole system under 
investigation here. As dynamical instabilities are signaled by the existence of ABMs 
with Im[u)] > 0, the complete absence of ABMs confirms the numerical observation 
of Sec. 3 that white holes are dynamically stable objects. Of course, the linearized 
Bogoliubov theory that is used in this Section is not able to capture the stationary 
Bogoliubov-Cerenkov density modulation pattern discussed in Sec. 3. 2, that appears as 
a consequence of nonlinear back-action effects. 

In a forthcoming paper this theory will be extended to more complex configuration 
involving a pair of black and white hole horizons: in agreement with [21, 22], the black 
hole laser phenomenon shows up as a series of discrete eigenmodes with a positive 
imaginary part of the frequency. Each of them corresponds to a zero of the determinant 
at the corresponding complex frequency Co. In a figure analogous to Fig. 11, such modes 
would appear as narrow blue spots corresponding to sharp minima of the logarithm. 

5.2. The S matrix and the white hole radiation 

In the previous subsection we have shown that no complex frequency eigenmodes appear 
in a white hole configuration. All the physics of quantum fluctuations can then be 
studied in terms of the input-output formalism of quantum optics [33, 34]: the central 
object of this approach is the S(u) matrix connecting the amplitude operators for the 
out-going Bogoliubov modes to those of the in-going modes at the same real frequency 
u. This is a direct extension to the white hole configuration of the Bogoliubov theory 
of acoustic black holes first introduced in [35] and fully developed in [18, 19, 25]. 

A sketch of the in-going and out-going modes for the white hole configuration 
was shown in Fig.2(a,b): in the low-energy sector below the cut-off frequency w max 
(defined as the maximum of the dispersion of the 112 mode), there are three in-going 
modes, labelled d m , u 1 ™ and u™, and three out-going modes, labelled d out , u° ut and n™ 1 . 
As it was discussed in previous works, the d and U\ out-going (in-going) modes have 
a positive Bogoliubov norm and enter in the input-output formalism as annihilation 
operators b dtUl (u) (d d ^ Ul (uj)). On the other hand, the negative norm u 2 modes enter the 
formalism as creation operators bl (oj) (0^ (u)). The input-output relations have the 

It can be shown that the number of such solutions is always two. More details on this issue as well as 
on the relation between this boundary condition and the in-going vs. out-going character of the modes 
will be given in future publications. 
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form: 



/ b d {u) 

KM ! = *( ■) 



( a d (u) 



(14) 



V 4 2 H 

The white hole radiation consists of a finite occupation of the out-going modes even for a 
vacuum state of the in-going ones. In this framework, it appears as a direct consequence 
of the mixing of creation and destruction operators by the S'-matrix (14). As usual in 
quantum optics [33, 34], such processes lead to a squeezed vacuum state as the output 
state of the modes after having scattered on the horizon. Exception made for some 
remarks at the end of Sec. 5. 4, our investigation is focussed on the case of a T = 
initial state where the in-going modes are assumed to be in their vacuum state defined 
by ad, Ul ,u 2 l vac ) — 0- The population of the out-going d and u\ modes is proportional to 
|5du 2 | 2 and (S'uiU.J 2 , respectively. Thanks to the unitarity of the S matrix, the population 
of the negative frequency partners is proportional to I^^J 2 = IS^] 2 + l^uj 2 , which 
confirms that phonons are created in pairs which carry no net energy. 

A simple way to estimate the matrix elements of S(u) is to look at the time evolution 
of incident wavepackets as done in Sec. 3, and to extract the reflection and transmission 
amplitudes from the intensity and phase of the out-going wavepackets. Alternatively, one 
can obtain the Bogoliubov modes at a given u> by solving the corresponding Bogoliubov- 
de Gennes differential equations: a numerical solution of this problem was performed 
in [18] for smooth configurations of the flow velocity and the speed of sound. Closed 
analytical formulas for the transmission and reflection amplitudes of a black hole in the 
opposite limit of a step-like horizon were obtained in [25] and then in [36]: imposing 
appropriate boundary conditions at the horizon x = 0, one is able to match the 
Bogoliubov plane-wave solutions within each homogeneous region for x ^ and then 
build a global solutions of the Bogoliubov problem. 

Here this last strategy is applied to a white hole configuration to obtain closed 
expressions for the matrix elements of S{u) in the low-frequency limit within the step- 
like horizon approximation. In particular, we shall focus our attention on the matrix 
elements for the negative frequency u 2 in-going mode that are involved in the white hole 
radiation: 

|C 132 I «y 12 ^ {^-C d ){vl-C 2 u f' 2 {c d + C u ) 1 

2c u (c d + v ){c u - c d ) w 
IW (16) 

c u [v + c d y 

This behavior of the low-energy S matrix elements is common to all in-going modes: all 
the matrix elements for out-going u\ and u 2 modes diverge as l/yju, while the ones for 
out-going d mode tends to a constant value. As a consequence, the spectral distribution 
of the quanta that are emitted into the modes into the white hole has the usual 1/u 
thermal form of Hawking emission, while the radiation leaving the horizon has a flat 
spectrum. This remarkable fact was first noticed in [18] as a result of a full Bogoliubov 
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calculation with a smooth horizon, and constitutes a crucial difference from the black 
hole case where the Hawking emission outside the horizon has a thermal 1/oj spectrum. 

For higher frequencies (but still below w max ), the S UlU2 matrix element quickly 
decreases to small values, while the S U2U2 one keeps a sizable value as it ends up 
describing simple reflection processes. For a smooth horizon of thickness a x ^> £, these 
coefficients are expected to tend to respectively and 1 with an exponential law of the 
canonical Hawking form 

\S^l\ = [l-eM-^M]- 1/2 (17) 
eM-"/2"H)\S U2U2 \, (18) 

where the specific value of the cut-off ojh = ^bTuI^ is determined by the analog of the 
surface gravity of the horizon configuration under examination [10]. 

Before proceeding with the calculation of the density correlations, it is interesting 
to point out a remarkable relation connecting the ^-matrix of black hole and white hole 
configurations related by a time-reversal symmetry [18, 19]: 

[S WH ]*= [S BH }-\ (19) 

with the convention that the d, ui, «2 modes of the white hole correspond after time- 
reversal (that exchanges the upstream and downstream regions) to the u, di, d<i ones of 
the black hole as defined in [25]. Making use of the 77-unitarity of S, that is S^rjS = rj 
with r] = diagfl, 1, —1], this leads to 

s BH = v [s WH ] T v . (20) 

5.3. Density correlations between the internal and the external regions 

The input-output approach reviewed in the previous section is the starting point for a 
calculation of the correlation function of density fluctuations. Such a calculation was 
first performed in [25] for the black hole case and is briefly reviewed in Sec. 5. 5. The 
results agreed with the numerical observations of [15]. Here we shall extend this theory 
to the case of white holes. 

Following the same lines as in [25], one can show that the density correlation 
between the internal and the external regions of the white hole far away from the 
horizon can be written in the form 

PCX) 

G-(2) ( x> a/) = l+ dcu 5Gf (u, x, x') (21) 
Jo 



with 

T7,n fiO^(u. X. X,') = 

4vr 



n SG '(u,x,x') = — \ \ j= — x 

Zl 7T ' * I /l rait rai+ I 



, I. /UiOut „,out 

l=u u u 2 \J\ U g,d u g,l 

x jl*S*W-*r*i»>)*] S* lU2 (co) S dU2 (co) + c.c.}. (22) 

Here, the point x is taken in the upstream region inside the white hole, while the point 
x' is taken in the downstream region outside the white hole. Of course, d~G (u>,x,x') 
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is non-zero only within the u < u max region where the u™ mode exists. As already 
mentioned, we are here restricting our attention to the zero temperature case, T = 0. 
Formula (22) clearly indicates that the process that is responsible for the correlation 
signal can be interpreted as the conversion of quantum fluctuations from the mode 
into out-going pairs of d ont and u° nt or d out and u^ ut real Bogoliubov excitations: the 
d ont quantum is emitted away from the white hole in the rightwards direction, while the 
u° U 2 quanta propagate upstream into the white hole. 

We have defined k^ 1U2 (u) as the value of the wavevector of the d,u\,U2 outgoing 
modes at the frequency u. The corresponding group velocities in the laboratory 
frame are defined as w° u ^ Ui n. The density components R k out ( w ) of the Bogoliubov 
modes are written in terms of the Bogoliubov Vk coefficients as Rk = £4 + Vk = 
[Hk 2 /2mQ(k)] 1 / 2 [9, 1]. In particular, note how R k out(u) goes as yfk^Huj) ~ y/cj 
in the low-frequency limit, while R^out^{uS) tend to finite values simply because both 

(u — y 0) = ±hz- this simple fact will play a crucial role in the following. 




Figure 12. Upper (a) and (b) panels, separated plots of the contribution of the d-U\ 
and d-U2 pairs to the correlation function of density fluctuations. Lower (c) plot, full 
analytical prediction (25) (red solid line). Numerical prediction evaluated at a late time 
fidt = 180 and rescaled to compensate for the spatial smoothening procedure (green 
dashed line). All the curves are cuts of the correlation function along the x' — 63.5^ 
line (magenta line in Fig. 7). The vertical dotted line indicates the position of the 
intersection of the cut line with the destructive interference line (26). 

Approximating the 5-matrix products with their low-energy form for a step-like 
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horizon configuration 
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and summing up the two terms, one is led to the final expression: 



(23) 
(24) 
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(25) 



Vz 



U 9,d 



X 



where = v + c d is the speed of low-energy out-going d phonons, while the low- 
frequency Mi 5 2 modes propagate at the zero-mode speed 1;°^ 2 ~ v z- 

The carrier wavevector of the fringe system is determined by the finite wavevector 
k u ont ~ kz of the out-going modes. In addition, its envelope shows a slower 
modulation at a wavevector u max /vz set by the ultra-violet scale ou max . The central 
minimum of the feature lies along the straight line 

r' 

—— = (26) 
v + c d 

and is due to a destructive interference between the u 1)2 contributions. This peculiar 
physics is illustrated in the cuts of the analytical pattern (25) that are shown in Fig. 12: 
panels (a,b) show the separate contributions of the Ui and u 2 modes: the carrier 
wavevector is set by kz, while the slow modulation is fixed by w max - The lower panel 
(c) shows their sum: destructive interference is responsible for the additional minimum 
in the envelope along the line (26). Given the number of approximations involved in 
the analytical calculations leding to (25), the qualitative agreement with the numerical 
results can be considered as quite good, in particular for what concerns the presence of 
the modulated envelope. 



5.4- Density correlations in the internal region: the growing checkerboard pattern 

The same framework can be applied to the correlations between points located in the 
internal region of the white hole far away from the horizon. In this case, auto-correlations 
in the U\ and u 2 outgoing modes are important, as well as cross-correlations between 
the two modes. Hence, the density correlation function is given by three contributions 



no 5G^\u, x, x') = 
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+ (xOx'). (27) 

The expressions for the low-frequency limit of the ^-matrix elements within the step-like 
horizon approximation are 
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(29) 



2c u (c (i + w )(c u -c (i )a; 
Inserting these expressions in (21), the integral over uj can be trivially performed. For 
the contribution of the u° ut — u™* cross-correlations and the sum of the m™' and u™ 1 
self-correlations, one obtains respectively: 



n 8G(?\x, x') 



{n-c 2 u ){c d -v ) 



2irc u (c d + v )(c u - c d ) 2 
x |(2w 2 - c\ - c 2 d ) cos[k z {x + x')] + 

+ 2^(c 2 -v 2 )(v 2 -c 2 u )sm{k z (x + x')]}l e , (30) 

noSGPM ^ + ^ oos[ib.(x-xO]/«. (31) 

u 1+ u 2 2irc u (c d + v ){c u - c d ) 

The trigonometric factors in (30) and (31) are responsible for the checkerboard pattern 

that was observed in Fig. 7: the former gives fringes parallel to the anti-diagonal 

x + x' = 0, while the former gives the fringes parallel to the main diagonal x — x' . 

The most interesting factor is however the integral 



— COS 



(32) 



vz 

which diverges when the infrared cut-off e — > 0. 

A way of coping with this divergence is to assume a scaling e oc 1/r for the infrared 
cut-off, r being the time elapsed since the formation of the horizon. In this way, I e 
becomes time-dependent: in particular, its value for x — x' follows a logarithmic growth 
with r. This fact is in agreement with the growth of the checkerboard pattern amplitude 
with time that was observed in the numerical simulations. In Fig. 9(a), we have tried to 
fit the numerical data with a logarithmic law (cyan line): the agreement is indeed quite 
good. 

As a further check of our interpretation, we have repeated the same calculation 
for the finite temperature T > case, which introduces an additional 1/uj factor into 
(32). The same strategy of imposing a infrared cut-off at 1/r then correctly explains the 
linear growth of the checkerboard pattern in time that was observed in the numerical 
calculation of Fig. 9(b). All these features will be discussed thoroughly in [32]. 
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5. 5. Analytic form of correlations in a black hole configuration 

For the sake of completeness and to better understand the peculiarities of the white 
hole emission, it is useful to briefly review the analytical form of the different features 
that appear in the correlation pattern of black holes [25] . The main ingredient of the 
calculation is again the S matrix that connects the amplitude operators for the out- 
going -u out , <i° ut , c?2 Ut Bogoliubov modes to those of the in-going u m , d™, d 2 n ones. From 
the relation (20), it is immediately clear that for all / = {u, d±, d 2 }, Si dl and S\ d2 scale 
as 1/y/ui, while Si u tends to a finite value in the low frequency limit. 

5.5.1. Density correlation between the internal and external regions The density 
correlation between the inside and outside regions of the black hole have the form: 

,^(2)/ /x 1 xr^ r-ftfcg ut ( a; )- R fc? ut ( a; ) 
n 5Gy(uj,x,x) = — > <^ x 

4 7T ^7^, L / 1 7, out 7,out| 
l=d 1: d 2 A/ \ V g,u U g,l I 

x ^H-W s^H S ud2 (u) + c.c.}, (33) 

Here, the point x > is located in the upstream region outside the black hole, white 
the point x 1 < is located in the downstream region inside the black hole. 

As in the white hole case, expression (33) clearly indicates that the process 
responsible for the correlation signal consists of quantum fluctuations incident on 
the horizon from the d 2 mode, which scatter into a correlated pair of real out- 
going Bogoliubov excitations either into the u and d 1 or into u and d 2 modes. As 
already mentioned, in the black hole configuration the wavevector kf nt of all out-going 
/ = {u, di,d 2 } modes tends linearly to zero in the zero frequency limit, which implies 
that all R™\u) oc y/w. 

For this reason, none of the black hole correlations shows any modulation and all 
features have a sharp tongue-like shape; within the low-frequency approximation of the 
S'-matrix mentioned above, one indeed obtains 

m ^ A^uwsinc (uw(^ - 

U V ' ) a out 7 .out /7T7T- V I 

l=d u d 2 ^\ V 9,u V Q,l W C uCd 

where sinc(a;) = sinx/x and the A[ coefficient has the following explicit form 

, [c^ (c u + vq) ( vl c 2 d \ 3/2 c u 

A dl[d2) = x — — —I—-— ——r^r (35) 



c d (c u - v ) \cl c\) c d + {-)c u 
in terms of the microscopic parameters of the system, v °^ = v + c u > is the group 
velocity of the w out mode escaping from the black hole. The group velocities v 2 < of 
the d° U 2 outgoing modes are very different from each other, 2=^0^ c d . This simple 
fact has the important consequence that the corresponding u-d^ and u-d1 ut features 
(iii) and (v) in the correlation pattern of Fig. 10(b) are spatially separated and do not 
interfere, as it was instead the case for a white hole case of Sec. 5. 3. Feature (iii) may 
be considered the main signature on the density correlation function of the Hawking 
emission from acoustic black holes. 
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5.5.2. Density correlations in the internal regions The correlation function for points 
both located inside the black hole, x, x' < contains a single contribution that originates 
from the di,d,2 modes, namely 



4:71 



„.out „.OUt I 
U g,di U 9,d2 I 



x S* dld2 (u)S d2d2 (uj)+c.c.) + (x^x'). (36) 

Within the same low-frequency expansion approximation, one obtains a tongue-shaped 
contribution of the form 



BQ cosine - d^)) 



with 



n 5Gi 2) (x, a;') ~ — y —^± + (x ^ x') . (37) 

2c d (c u -v ) \c 2 u c 2 u J 

that correctly reproduces feature (iv) of Fig. 10(b). In particular, differently from the 
white hole case, no any infra-red divergence appears in the integral, which confirms the 
stationarity in time of the Hawking radiation pattern. 



In summary, in this Section we have developed a Bogoliubov theory of quantum 
fluctuations in a white hole configuration within the step-like approximation. This 
theory provides analytical formulas for the physical observables that were considered in 
the numerical simulations reported in the previous Section. In contrast to the thermal 
character of the Hawking radiation from a black hole, the spectral distribution of the 
emitted radiation outside the white hole turns out to be flat. The steady growth of 
the density fluctuations in time is signalled by infrared divergences in the Bogoliubov 
calculation. By imposing a suitable cut-off to the integrals one is able to reproduce 
the numerically observed trend. Inclusion of the terms describing the back-action 
of fluctuations onto the condensate is presently in progress and is expected to shine 
light on the effect of quantum fluctuations on the stability of the acoustic white hole 
configuration. 



6. Conclusions 



In this paper we have reported a comprehensive investigation of the physics of acoustic 
white hole configurations. We have focussed our attention on a simplest model that 
is amenable to analytical and numerical work. Using the Gross-Pitaevskii equation for 
the condensate dynamics, we have established the dynamical stability of the white hole 
at the mean-field level: no exponentially growing perturbation of the horizon appears 
when this is hit by an incident wavepacket of Bogoliubov phonons. A nonlinear back- 
action effect due to the perturbation of the horizon by the incident wavepacket is visible 
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as a stationary pattern of Bogoliubov-Cerenkov waves by the perturbed horizon; the 
amplitude of this pattern scales as the square of the incident wavepacket one. 

Quantum fluctuation are included within a truncated Wigner numerical approach; 
the physics of an infinitely long condensate is reproduced by imposing absorbing 
boundary conditions to fluctuations. Signatures of the white hole emission by the 
horizon are identified in the correlation function of density fluctuations. Differently 
from the Hawking emission from black holes, the white hole emission does not consist 
of low-energy excitations only, but involves high wavevector modes. As a result, the 
density correlation pattern shows fastly oscillating features on the scale of the healing 
length. Furthermore, some of these features are not stationary in time, but indefinitely 
grow with time with a logarithmic (linear) law for a zero (finite) initial temperature. 
All the features found in the numerical observations are fully recovered and physically 
interpreted by an analytical calculation within a step-like approximation for the horizon. 
Interestingly, this theory predicts a flat spectral distribution of the radiation emitted 
outside the acoustic white hole: this is in strike contrast to the thermal character of 
Hawking radiation from black holes. 

Besides its intrinsic interest from the point of view of quantum fluctuations in 
condensed matter systems, our results may open important avenues in gravitational 
physics and quantum field theory in curved space-times: in spite of being related to 
each other by a simple time reversal symmetry, the physical properties of the quantum 
emission from black and white holes turn out to be significantly different, starting from 
the thermal vs. flat frequency distribution of the emission. On one hand, Hawking 
radiation from black holes is concentrated in the low- A; hydrodynamical modes whose 
dispersion has a linear, relativistic-like behaviour; for this reason, Hawking radiation 
has turned out to be very robust against modifications of the ultraviolet behaviour of 
the dispersion relation. On the other hand, white hole radiation is intrinsecally not an 
hydrodynamical process: because of the large associated blue shift, it is in fact highly 
sensitive to the ultraviolet behaviour of the Bogoliubov dispersion relation. On this 
basis, we expect that white hole radiation may then provide an even better magnifying 
glass to test the intimate microscopic structure of the quantum vacuum. 

Acknowledgements 

We are indebted to N. Pavloff and M. Rinaldi for stimulating discussions. A. F. wishes 
to thank Generalitat Valenciana for financial support. A.F. and CM. are supported in 
part by MICINN grant FIS2008-06078-C03-02. CM. is indebted to the Spanish MEC 
for a FPU fellowship. CM. would like to thank the Universita di Bologna and the BEC 
Center in Trento for their charming hospitality during the elaboration of this work. 

References 

[1] L. Pitaevskii and S. Stringari, Bose-Einstein condensation (Clarendon Press, Oxford, 2003). 



Acoustic white holes in atomic BECs 



31 



[2] I. Carusotto, S. X. Hu, L. A. Collins, A. Smerzi, Phys. Rev. Lett. 97, 260403 (2007). 
[3] P. Engels and C. Atherton, Phys. Rev. Lett. 99, 160405 (2007). 

[4] L. Fallani, L. De Sarlo, J. E. Lye, M. Modugno, R. Saers, C. Fort, and M. Inguscio, Phys. Rev. 
Lett. 93, 140406 (2004). 

[5] B. Wu and Q. Niu, Phys. Rev. A 64, 061603(R) (2001); A. Smerzi, A. Trombettoni, P. G. 

Kevrekidis, and A. R. Bishop, Phys. Rev. Lett. 89, 170402 (2002); M. Modugno, C. Tozzo, 

and F. Dalfovo, Phys. Rev. A 70, 043625 (2004). 
[6] V. Hakim, Phys. Rev. E 55, 2835 (1997); N. Pavloff, Phys. Rev. A 66, 013610 (2002). 
[7] T. Frisch, Y. Pomcau, S. Rica. Phys. Rev. Lett. 69 1644 (1992); T. Winiecki, B. Jackson, J. F. 

McCann and C. S. Adams, J. Phys. B: At. Mol. Opt. Phys. 33, 4069 (2000); A. Kamchatnov 

and L Pitacvskii. Phys. Rev. Lett. 100 160402 (2008). 
[8] A. M. Leszczyszyn, G. A. El, Yu. G. Gladush, and A. M. Kamchatnov, Phys. Rev. A 79, 063608 

(2009). 

[9] Y. Castin, in 'Coherent atomic matter waves', Lecture Notes of Les Houchcs Summer School, 
p. 1-136, edited by R. Kaiser, C. Westbrook, and F. David, EDP Sciences and Springer- Verlag 
(2001) 

[10] W.G. Unruh, Phys. Rev. Lett. 46, 1351 (1981). 

[11] Artificial Black Holes, edited by M. Novello, M. Visser, and G. Volovik (World Scientific, River 
Edge, 2002). 

[12] N. D. Birrell and P. C. W. Davies, Quantum fields in curved space (Cambridge University Press, 
Cambridge, UK, 1982). 

[13] S. W. Hawking, Nature (London) 248, 30 (1974); S. W. Hawking, Commun. Math. Phys. 43, 199 
(1975). 

[14] R. Balbinot, A. Fabbri, S. Fagnocchi, A. Recati, and I. Carusotto, Phys. Rev. A 78, 021603 (2008). 

[15] I. Carusotto, S. Fagnocchi, A. Recati, R. Balbinot, and A. Fabbri, New J. Phys. 10, 103001 (2008). 

[16] O. Lahav, A. Itah, A. Blumkin, C. Gordon, J. Steinhauer, preprint arXiv:0906.1337. 

[17] U. Lconhardt, T. Kiss, and P. Ohbcrg, Phys. Rev. A 67, 033602 (2003). 

[18] J. Machcr and R. Parentani, Phys. Rev. A 80, 043601 (2009). 

[19] J. Machcr and R. Parentani, Phys. Rev. D 79, 124008 (2009). 

[20] C. Barcelo, A. Cano, L. J. Garay, and G. Jannes, Phys. Rev. D 74, 024008 (2006). 

[21] A. Coutant and R. Parentani, Phys. Rev. D 81, 084042 (2010). 

[22] S. Finazzi and R. Parentani, preprint arXiv: 1005.4024. 

[23] L. J. Garay, J. R. Anglin, J. I. Cirac, and P. Zollcr, Phys. Rev. A — bf 63, 023611 (2001); P. Jain, 

A. S. Bradley, and C. W. Gardiner, Phys. Rev. A 76, 023617 (2007). 
[24] S. Corley and T. Jacobson, Phys. Rev. D 59, 124011 (1999). 
[25] A. Recati, N. Pavloff, I. Carusotto, Phys. Rev. A 80, 043603 (2009). 
[26] R. Balbinot, A. Fabbri, S. Fagnocchi, and R. Parentani, Riv. Nuovo Cim. 28, 1 (2005). 
[27] C. Barcelo, S. Liberati, and M. Visser, Living Rev. Relativity 8, 12 (2005). 

[28] S. Weinfurtner, private communication. For details on the experimental set-up, see S. Weinfurtner, 
E. W. Tedford, M. C. J. Prentice, W. G. Unruh, G. A. Lawrence, preprint arXiv:1008.1911. 

[29] G. Rousscaux, C. Mathis, P. Mai'ssa, T. G. Philbin and U. Leonhardt, New J. Phys. 10, 053015 
(2008). 

[30] A. Sinatra, C. Lobo, Y. Castin, J. Phys. B 35, 3599 (2002). 
[31] I. Carusotto, R. Balbinot, A. Fabbri, A. Recati, Eur. Phys. J. D 56, 391 (2010). 
[32] C. Mayoral Saenz, PhD thesis at the Universidad de Valencia, expected for 2011. 
[33] D.F. Walls and G.J. Milburn, Quantum Optics (Springer- verlag, Berlin, 1994). 
[34] C. Ciuti and I. Carusotto, Phys. Rev. A 74, 033811 (2006). 

[35] U. Leonhardt, T. Kiss, and P. Ohbcrg, J. Opt. B: Quantum Semiclassical. Opt. 5, S42 (2003). 
[36] A. Fabbri and C. Mayoral, arXiv: 1004.4876; C. Mayoral, A. Fabbri and M. Rinaldi, 
arXiv:1008.2125. 



